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i In this talk I discuss the theoretical issues associated with lattice calculations of isospin 

breaking corrections to hadronic matrix elements. I concentrate on the calculation of QCD 
isospin breaking effects for the K£2 and K£3 decay rates and illustrate the recent lattice 
results obtained by the RM123 collaboration. 
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1 Introduction 

oo 

OO Isospin is an approximate symmetry of fundamental interactions. The two lightest 

quarks, the up and the down, have different masses and different electric charges 
but it happens that their mass difference is a small parameter with respect to the 

fTj confinement scale Kqcd and that electromagnetic isospin breaking effects are also 

small on hadronic observables because at low energy a em is much smaller than a s . 
For these reasons most of the theoretical predictions on phenomenologically relevant 



hadronic observables have been derived by assuming the exact validity of isospin 
symmetry. This is also the case for most of the non-perturbative theoretical predic- 
tions on hadronic matrix elements obtained over the years by performing lattice QCD 
simulations. 

In this talk I will concentrate on the lattice calculation of isospin breaking correc- 
tions to the hadronic matrix elements entering the leptonic (K£2) and semileptonic 
(K£3) decay rates of the kaons. These are among the golden-plated quantities for 
lattice QCD calculations: these involve a single hadron into initial and/or final states. 
Furthermore, two and three-point functions of pseudoscalar mesons interpolating op- 
erators can be obtained with high statistical accuracy by using a reasonably small 
number of gauge configurations. 
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In order to highlight the importance of a non-perturbative first principle calcu- 
lation of the isospin breaking effects on these quantities I quote below the FLAG PQ 
averages of the lattice results for the ratio of the kaons (F K ) and pions (F^) decay 
constants and for the form factor parametrizing the semileptonic decay of a kaon into 
apion(,/f% 2 )), 



FLAG 



= 1.193(5) 



/r(o) 



FLAG 



= 0.956(8) . 



(1) 



These results assume exact isospin symmetry. One may wonder whether the quoted 
uncertainties (0.5% and 0.8% respectively) can be considered reliable in view of all 
the systematics sources of uncertainties that have to be taken into account when a 
lattice calculation is performed (unquenching of the heavy flavours, chiral extrapola- 
tions, finite volume effects, renormalization, continuum extrapolations, etc.). For a 
detailed discussion of all the systematics affecting lattice calculations of K£2 and K£3 
decay rates see, for example, the contributions of J. Laiho and of A. Jiittner to these 
proceedings. Here I argue that the figures quoted above are in fact reliable, by using 
the argument of the FlaviaNet working group [2] (see also [3] and the contribution 
of V. Cirigliano to these proceedings) based on the first-row CKM unitarity test. 
Indeed, by neglecting V u b because it is smaller than the uncertainties on the other 
quantities, by using the experimental determinations of the pions and kaons leptonic 
decay rates together with the experimental determination of the K£3 decay rate and 
by using V u a as extracted by the combined experimental and theoretical analysis of 
^ 20 super-allowed nuclear /3-decays, 



V US F K 



VudFn 



exp 



0.2758(5) 



Kjf^O) ~" = 0.2163(5) , 



exp 



IK 



ud | 



\exp 



0.97425(22) , 



\V u d\ 2 + \Vusf 



one is able to get a "measurement" of the hadronic unknowns, 



K 



CKM 



1.1919(57) 



ff*(P) 



CKM 



0.9595(46) . 



(2) 



(3) 



The numbers above have been obtained without using lattice inputs and have been 
corrected for isospin breaking effects by using the effective field theory calculations of 
refs. [HE]. Concerning the QCD isospin breaking effects, i.e. those coming from the 
difference of the up and down quark masses in absence of electromagnetic interactions, 
the authors of refs. [H E] obtain 



Fk+/F^ 
, F K /F V 



x QCD, X pt 

1 = -0.0022(6) 



7r*V) 



v QCD, X pt 

1 I = 0.029(6) . (4) 
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There are some considerations that follow from the comparison of eqs. ([I]) with 
eqs. ([3]) and eqs. The uncertainties quoted by lattice calculations are in fact 

reliable and, at least for these gold-plated quantities, lattice QCD calculations entered 
the 1% level of accuracy. Furthermore, isospin breaking effects on these quantities 
are as big as the current lattice errors and cannot be neglected in future lattice 
calculations, it would be useless otherwise to put efforts in further reducing lattice 
errors on the isosymmetric quantities. Isospin breaking effects are particularly large in 
the case of the K£3 decay rate. By assuming SU (3)-flavour symmetry the form factor 
at vanishing momentum transfer is constrained to be equal to one by the conservation 
of the vector current and isospin breaking effects are of the same order of magnitude 
of the deviation of /^(O) from this normalization condition. 

In the rest of this talk I describe the method recently developed by the RM123 
collaboration [6] in order to calculate on the lattice QCD isospin breaking effects and 
I discuss the results obtained by applying this method. In the last part of the talk 
I also briefly comment on the lattice calculation of electromagnetic isospin breaking 
effects. 



2 Leading QCD isospin breaking effects 

The RM123 collaboration has recently proposed [6] a method to calculate leading 
QCD isospin breaking effects on the lattice by starting from gauge configurations 
generated within the isosymmetric QCD theory (m^ = m u ). The method is based on 
the perturbative expansion of the path-integral with respect to the small parameter 
Am ud = (m d - m u )/2, 

{0) „ JD0O(1 + Am ud S) e~ s ° = (O)o + Am ud (OS) 
I D<t> (1 + Am ud S) e- s o 1 + Am ud (S) 

= (0)o + Am ud (OS) , (5) 

where (-) represent the vacuum expectation value in the isospin symmetric theory and 
S = J2 X \y* u ~ dd](x) is the isospin breaking term. The correction in the denominator 
vanishes, (5% = 0, because of isospin symmetry. In order to take into account first 
order insertions of S into lattice correlators, one can start by the given observable 
into the full isospin broken theory, integrate out the quark fields (take the fermionic 
Wick contractions) and expand the quark propagators at fixed QCD background field 
with respect to Am ud . Diagrammatically this amounts to 

u 

= — > + Amy — (x) 1 , 

d 

= > Am ud — ® 1 , (6) 
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Figure 1: Left panel: combined chiral and continuum extrapolations of AM|/Am U( |. Right 
panel: combined chiral and continuum extrapolations of 5 Fx / Amy. Black points correspond to 
the coarser lattice spacing, a = 0.098 fm, dark magenta points correspond to a = 0.085 fm, green 
points to a = 0.067 fm and blue points to a = 0.054 fm. Red lines are the results of the continuum 
extrapolations. 



where the propagators with an insertion are simply obtained by calculating the square 
of the propagator in the isosymmetric theory. By using this simple recipe (see ref. [6] 
for further details and for all the unexplained notation) is easy to show that pion 
masses and decay constants do not get a first-order QCD isospin breaking correction. 
Concerning pseudoscalar-pseudoscalar two-point correlators of the kaons, one gets 

C k +k- (t) = - <^ = - ^> - Am ud <^> + 0(Am ud ) 2 , 



m ud ) 



(7) 



where the strange quark propagator has been drawn in red and where the same color 
code of eqs. ^ has been used for the up and down propagators. By studying the 
euclidean time behavior of the ratio of correlators appearing on the right-hand sides 
of the previous equations, 



AC KK (t) 
Amy 



(8) 



it is possible to calculate the leading QCD isospin breaking effects on kaon masses and 
decay constants. In Figure[I]l show the combined chiral and continuum extrapolations 
of the lattice data of ref. [6j from which the following results 



Amy 



QCD 



(MS, 2GeV) 



M 2 K0 



M 2 K+ 



m d 



m,, 



QCD 



= 2.57(8) x 10 3 MeV 



4 



K 



Am, 



QCD 



(MS, 2GeV) 



F K — F K + 
F K (m d - m u ) 



QCD 



3.3(3) x 10" 3 MeV~\(9) 



are obtained. The results above have been obtained by including in the generation 
of the gauge configurations the determinants of the two lightest quarks, i.e. in the 
so-called nj — 2 theory. By taking the ratio of the previous two results it is possible 
to express the correction to the K£2 decay rate, and consequently to F K fF^ (given 
the fact that SF^ = at this order), as a function of the difference of the neutral and 
charged kaon masses induced by the quark mass difference in QCD. 

It is important to realize that the measured mass difference of the kaons is due 
to QED as well as to QCD isospin breaking effects and that the first are expected 
to be not negligible with respect to the latter. Furthermore, the separation of QED 
from QCD isospin breaking effects is arbitrary and prescription dependent because of 
the renormalization of the quark masses and of the strong coupling constant required 
in order to absorb the divergences generated by electromagnetic interactions. QED 
isospin breaking effects on hadronic quantities can also be calculated on the lattice. 
The first pioneering calculation has been performed in ref. [7J where lattice QED 
has been treated in its "non-compact" formulation. Later on, several lattice collab- 
orations provided calculations of the electromagnetic mass splittings of the lightest 
hadrons, see for example refs. JH [9J UHl E] for recent works on the subject. 

In ref. [12] the RM123 collaboration is going to provide a method to calculate 
QED+QCD isospin breaking effects by starting from isosymmetric QCD gauge con- 
figurations while, for the time being, they have used the following prescription 



M 2 



K° 



Ml + 



QCD 



M 2 



M 2 



exp 



(l + e 7 ) M 2 



M 2 



exp 



6.05(63) x 10 3 MeV 2 



0.7(5) 



(10) 



with e 1 taken from [lj in order to define the QCD part of the kaons mass difference. 
The physical meaning of the previous prescription comes from the observation that, 
according to the well known Dashen's theorem, the mass splittings of the kaons and of 
the pions have to be equal into the (m 8 = ma = m u = 0) chiral limit. The (relatively) 
small parameter e 1 parametrizes the violation of the Dashen's theorem and can be 
calculated in a given renormalization scheme by using lattice results. By using this 
prescription, and the results of eqs. ([9]), the RM123 collaboration quotes 



F K +/F^ 
F K /F n 



QCD 



-0.0039(3)(2) 



[M 2 



M 2 K+ 



\QCD 



6.05 x 10 3 MeV 2 



The result above is of the same order of magnitude, though higher, of the chiral 
perturbation theory estimate of ref. [5] (see eqs. ^ above) and strongly depend upon 
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Figure 2: Left panel: results for f^{q 2 ). Right panel: results for Sfff 7r (q 2 ). The data are 
obtained at fixed lattice spacing a = 0.085 fm and at fixed am^ = 0.0064. 



the value of e 1 . It is for this reason that is particularly important to perform combined 
lattice calculation of QCD+QED isospin breaking effects as it will be done in ref. [12]. 

In the last part of this talk I briefly discuss the issues related to the calculation 
of QCD isospin breaking corrections to the K£3 decay rate. By using the same 
techniques that lead to eq. ([8]), in ref. [6j it have been derived all the formulae needed 
in order to perform a lattice calculation of f+ +7T °(q 2 )/f+ 07T (q 2 ). In this case the 
deviations of the numerator and of the denominator from the common isosymmetric 
limit f K7T (q 2 ) are different. In particular, the formulae for f+ +7T (q 2 ) present quark 
disconnected (gluon connected) diagrams that are difficult to calculate on the lattice. 
These disconnected contractions are responsible of the tt — 77 mixing and are expected 
to give a large "resonant" contribution to the deviation of f+ +7T °(q 2 )/f+ 07T (q 2 ) from 
one. In ref. [6] these disconnected contributions have not been calculated while it has 
been provided a preliminary estimate of the relative deviation 

QCD 



/f-(0) 



I QCD 



0.85(18)(1) x 10" 



x 



[M 2 K0 - M 2 K+ 
6.05 x 10 3 MeV 2 



(12) 

In Figure [2] are shown the corresponding results for the isosymmetric form factor and 
for the correction as a function of the momentum transfer. The argument concerning 



the importance of disconnected diagrams is confirmed by the result of eq. (12) that is 



more than two orders of magnitude smaller of the chiral perturbation theory estimate 
of the physical correction performed in ref. [4] and quoted into eqs. ([3]) above. 

3 Outlooks 



In this talk I have argued that lattice QCD calculations are nowadays sufficiently 
accurate to be able to provide results for the K£2 and K£3 decay rates with sub- 
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percent overall uncertainties. Furthermore, the lattice errors on these quantities can 
be considered reliable, as emerges from the comparison of world average lattice results 
with the corresponding quantities obtained from the first-row CKM unitarity test. 

At this level of accuracy isospin breaking corrections cannot be neglected anymore 
because, at least for these quantities, these are of the same order of magnitude of the 
quoted uncertainties. Isospin breaking corrections can be estimated by recurring 
to chiral perturbation theory but can also be calculated from first principle lattice 
calculations. The RM123 collaboration has recently obtained encouraging results for 
the QCD isospin breaking corrections to the K£2 decay rate and preliminary results 
for the K£3 decay rate. 

QED isospin breaking corrections are expected to be as important as QCD ones 
and are required also for the calculation of the QCD effects in order to cope with the 
prescription dependent separation of electromagnetic and QCD contributions to the 
hadrons mass splittings. The QED corrections to the hadron spectrum can also be 
calculated on the lattice and several lattice collaborations have provided phenomeno- 
logically relevant results in the last years. 

Concerning QED corrections to weak matrix elements, I want to note in pass- 
ing that the problem is much more involved with respect to the case of the hadron 
spectrum. Already in the case of the QED corrections to the K£2 decay rate, elec- 
tromagnetic gauge invariance and cancellation of infrared divergences require the 
inclusion of correlators in which photons connect quark and lepton propagators. For 
this reason it is not possible to factorize the leptonic tensor in the calculation of the 
decay rate and more theoretical work is required before we shall be able to provide 
a first principle lattice calculation of the QED+QCD isospin breaking corrections to 
K£2 and K£3 processes. 

I warmly thank my colleagues of the RM123 collaboration for the enjoyable work 
of the last two years: the lattice results and the theoretical arguments discussed in 
this talk are the outcome of this work. I am grateful to the working group conveners 
and to the organizers of the CKM2012 conference for the opportunity of giving this 
talk. 
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